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V. Rotival,i'2'El K. Bennaceur,3^4,[t] ^nd T. Duguet^'^.s.ji] 

' DPTA/Service de Physique Nucleatre - CEA/DAM Ile-de-France - BP12 - 91680 Bruyeres-le-Chdtel, France 

''National Superconducting Cyclotron Laboratory, 1 Cyclotron Laboratory, East Lansing, MI 48824, USA 

^ Universite de Lyon, F-69003 Lyon, France; Institut de Physique Nucleaire de Lyon, 

CNRS/IN2P3, Universite Lyon 1, F-69622 Villeurbanne, France 

^ CEA, Centre de Saclay, IRFU /Service de Physique Nucleaire, F-91191 Cif-sur-Yvette, France 

''Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA 

(Dated: March 22, 2009) 

The analysis method proposed in Ref. [l| is applied to characterize halo properties in finite many- 
fermion systems. First, the versatility of the method is highlighted by applying it to light and 
medium-mass nuclei as well as to atom-positron and ion-positronium complexes. Second, the de- 
pendence of nuclear halo properties on the characteristics of the energy density functional used 
in self-consistent Hartree-Fock-Bogoliubov calculations is studied. We focus in particular on the 
influence of (i) the scheme used to regularize/renormalize the ultra-violet divergence of the local 
pairing functional, (ii) the angular-momentum cutofi' in the single-particle basis, as well as (iii) the 
isoscalar efi'ective mass, (iv) saturation density and (v) tensor terms characterizing the particle-hole 
part of the energy functional. It is found that (a) the low-density behavior of the pairing functional 
and the regularization/renormalization scheme must be chosen coherently and with care to provide 
meaningful predictions, (b) the impact of pairing correlations on halo properties is significant and 
is the result of two competing effects, (c) the detailed characteristics of the pairing functional has 
however only little importance, (d) halo properties depend significantly on any ingredient of the 
energy density functional that infiuences the location of single-particle levels; i.e. the effective mass, 
the tensor terms and the saturation density of nuclear matter. The latter dependencies give insights 
to how experimental data on medium-mass drip-line nuclei can be used in the distant future to 
constrain some characteristics of the nuclear energy density functional. Last but not least, large 
scale predictions of halos among all spherical even-even nuclei are performed using specific sets of 
particle-hole and particle-particle energy functionals. It is shown that halos in the ground state 
of medium-mass nuclei will only be found at the very limit of neutron stability and for a limited 
number of elements. 

PACS numbers: 2.10.Gv, 21.10.Pc, 21.30.Fc, 21.60.Jz 
Keywords: halo; EDF; HFB; pairing 
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I. INTRODUCTION 



The formation of halos is a quantum phenomenon 
caused by the possibility for non-classical systems to 
expand in the classically forbidden region j2i S S M, 
Q. Indeed, weakly-bound systems can extend well be- 
yond the classically-allowed region, as has been theoreti- 
cally predicted or experimentally observed for molecules 
(3He-3He-39K [7|, ''Hea BiEl, ^He^Hes [uL.), atom- 
positron complexes (e+Be, PsLi+, PsHe"'"...) [l^l and hy- 
pernuclei (jji) fl3- In nuclear physics, where the study 
of halos was initiated, efforts are still devoted to reach a 
better understanding of the structure and reaction prop- 
erties of such exotic systems. For instance, the existence 
of halos in borromean systems or excited states of mirror 
nuclei still raises questions [14|. In light nuclei, it was 
found that a cluster picture is at play, and known halo 



•Electronic address: 'vincent.rotival@polytechn ique.org| 
^Electronic address: bennaceur@ipnl.in2p3.fr 
■••Electronic address: |thonias.duguet@cea.fr| 



systems are accurately described by two- [il,[ia| or three- 
body [1^, ll3, Ha, llMl models, where one or two nucleons 
evolve around a tightly-bound core. This leads to a clas- 
sification in terms of one-nucleon halos (^^Be (201. l2ll |23| . 
19 C [2I [13, i^Ne Jp, d^...) and two-nucleon halos 
(^He [l|, "Li [13, [11^.). 

Since the discovery of the anomalous cross-section of 
11 Li [23, [2^1, one of the compelling questions relates to 
the existence of a mass limit beyond which the forma- 
tion of halos is inhibited. On the proton-rich side, it is 
believed that the Coulomb interaction prevents the for- 
mation of halos beyond Z « 10 [23]. However, this could 
be put into question as non-trivial effects may come into 
play [30[. From a theoretical standpoint, halos in medium 
to heavy mass nuclei can be studied through relativis- 
tic or non-relativistic Hartree-Fock-Bogoliubov calcula- 
tions [3ll [32] performed in the context of energy den- 
sity functional (EDF) methods. On the experimental 
side, the next generation of radioactive ion beam facil- 
ities (FAIR at GSI, RIBF at RIKEN, REX-ISOLDE at 
CERN, SPIRAL2 at GANIL...) might be able to as- 
sess the position of the neutron drip-line up to about 
Z « 26 ^3^. Although this would be an astonishing 



accomplishment, it will not allow the study of most of 
potential medium-mass halos. Still, the (distant) fu- 
ture confrontation of theoretical results with experimen- 
tal data will provide crucial information that can be used 
to constrain theoretical models. 

One difficulty resides in the absence of tools to char- 
acterize halo properties of finite many-fermion systems 
in a quantitative way. Light nuclei constitute an excep- 
tion considering that the quantification of halo proper- 
ties in terms of the dominance of a cluster configura- 
tion and of the probability of the weakly-bound clusters 
to extend beyond the classical turning point is well ac- 
knowledged [2^,[3l,[33,[3a| • Existing definitions and tools 
applicable to systems constituted of tens of fermions are 
too qualitative, the associated observables are incomplete 
and have led to misinterpretations in nuclei jl|. 

To improve on such a situation, a new quantitative 
and model-independent analysis method was proposed 
in Ref. [ll|. The method uses universal properties of 
the internal one-body density to extract, in a model- 
independent fashion, the part of the density that can be 
identified as a halo. Two criteria have been introduced 
to characterize halo systems in terms of (i) the average 
number of fermions participating in the halo and (ii) the 
infiuencc of the latter on the system extension. The re- 
sults deduced from EDF calculations of medium-mass 
nuclei have underlined the likely formation of a collec- 
tive halo at the neutron drip-line of chromium isotopes. 
The neutron density of those nuclei displays a spatially 
decorrelated region built out of an admixture oi £ = 
and £ — 2 overlap functions (orbitals). The significant 
contribution of orbitals with orbital angular-momentum 
£ = 2 is at variance with the standard picture in light 
nuclei [13, [s^. 

Given such an analysis method, important questions 
can now be addressed. The versatility of the method 
must be tested which we do in the first part of our study 
by applying it to many-fermion systems of different scales 
computed with various many-body techniques : atom- 
positron/ion-positronium complexes on the one hand and 
light/heavy nuclei on the other. In the second part of the 
article, we focus on the influence of (i) the scheme used to 
regularize/renormalize the ultra-violet divergence of the 
local pairing functional, (ii) the angular-momentum cut- 
off in the single-particle basis, as well as (iii) the isoscalar 
effective mass, (iv) saturation density of nuclear matter 
and (v) tensor terms [39] characterizing the particle-hole 
part of the energy functional. More generally, the dif- 
ferent ways pairing correlations impact halo nuclei are 
studied, e.g. the anti-halo effect [40, [4l| or the po- 
tential decorrelation of ^ = orbitals from the pairing 
field [Ulillil. 

The present paper is organized as follows. The anal- 
ysis method proposed in Ref. [l|] is brieffy recalled in 
Sec. In] whereas its versatility is highlighted in Sec. IIIBI 
In Sec. mil technical aspects of Skyrme-EDF calculations 
are provided and the dependence of halo predictions on 
some their ingredients is pointed out. Section llVI is de- 



voted to discussing the effect of pairing correlations on 
the formation of halos. Then, the sensibility of halo prop- 
erties to the characteristics of the particle- hole part of the 
functional is studied. A large scale study of potential ha- 
los among all spherical medium-mass nuclei is proposed 
in Sec. |Vl In conclusions, we discuss how halo systems 
could help constraining the nuclear EDF and to which ex- 
tent the present results can be related to data generated 
in the distant future by radioactive ion beam facilities. 



II. CHARACTERIZATION OF HALO SYSTEMS 
A. Analysis method 

Anew quantitative analysis-method of halos in finite 
many-fermion systems was proposed in Ref. [l|. The 
starting point is a model-independent decomposition of 
the internal one-body density prij (r) of spherical many- 
fermion systems in terms of spectroscopic amp litudes 
(/?y(r ) and their radial components ifn^e^j^ir) |45l.l46ll47| 
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The appearance of a halo in the A^-body system, i.e. 
the part of the density that is spatially decorrelated from 
an a priori unknown core, was shown in Ref. py] to be re- 
lated to the existence of three typical energy scales in the 
excitation spectrum of the (TV— l)-body system. From a 
practical viewpoint, and using the internal one-body den- 
sity as the only input, the method allows the extraction 
of the radius rg beyond which the halo, if it exists, is 
located. 

With the radius tq at hand, two quantitative halo 
factors are introduced. First, the average number of 
fermions participating in the halo can be extracted 
through 



+ 00 
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Second, the contribution of the halo region to the root 
mean square radius of the system is also extracted 
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B. Versatility of the method 

To illustrate the analysis method and its versatility, 
we now apply it to the results of many-body calculations 
performed for three different systems: light nuclei studied 
through coupled-channels calculations Iq, ll8[, medium- 
mass nuclei described through single-reference energy 





Core 




Total 




Mialo 5_Rhalo [fm] 


A^halo 


(J-Rhalo [fm] Rr.m.s. [fm] 


"Be 


0.000 0.000 
0.000 0.000 


0.66 X 10" 
0.270 


"3 0.74 X 10^3 2.487 
0.394 2.908 



Table 1: Values of A^haio and (5i?haio for ^^C (stable) and "Be 
(one-neutron halo). The criteria are applied to results ob- 
tained through coupled-channels calculations [l6l.ll8|]. 



density functional calculations and atom-positron/ion- 
positroniuni complexes computed with the fixed-core 
stochastic variational method [48|, |43, [50, IMl l52| • 



1. Light nuclei 

To check the consistency of the method in a situation 
where core and halo densities are explicitly computed, 
the values of A'haio and (5-Rhaio have been extracted from 
coupled-channels calculations of light nuclei. This will 
serve a basis for cases where the total density only is 
accessible. 

Calculations are performed for a core-|-neutron system 
for which the internal dynamics of the core is taken into 
account [la, uM and the total Hamiltonian reads as 



Htct — Hr, 



- Trcl + 1^1- 



(4) 



To provide adequate nuclear quadrupole couplings, a 
deformed Woods-Saxon potential in the core rest frame 
is considered 
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where j3 is the core quadrupole deformation. The total 
wave- function is expanded in a basis of eigenstates of the 
total angular momentum using a separation of the core 
internal motion, with eigenstates associated to the ener- 
gies ti, from the neutron relative motion. The resulting 
coupled-channels equations for the loosely bound nucleon 
wave function '^i read 



(T„i + e, - i;) *. + ^ y. 



^3 *J = 



(7) 



and are solved in a Sturmian basis. 

Calculations have been performed for two nuclei: the 
well-established one-neutron halo ^^Be, whose total den- 
sity and core+halo decomposition are plotted in Fig. [U 
and the stable nucleus ^^C as a control case. The results 
are summarized in Tab.|Tl where the two criteria are eval- 
uated for core and total densities. For ^^Be, ro is found 
to be compatible with the condition to have one order 
of magnitude difference between core and halo densities 
at that radius |1||. The tail-to-core ratio is slightly dif- 
ferent from ten, partly because the core is represented 
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Figure 1: (Color Online) Core, halo and total neutron de nsity 
of "Be obtained from coupled-channels calculations [la. [18]. 
The value of ro and the tolerance margin are indicted by ver- 
tical lines. 



by a gaussian profile with the wrong asymptotic. In any 
case, the ideal value of ro still lies within the allowed the- 
oretical error. This validates the method on a realistic 
system. The halo parameter A'haio shows that around 0.3 
neutron reside in the decorrelated region in average. The 
reason why one only finds a fraction of a neutron within 
the halo region is because the wave function of the " halo 
nucleon" lies partly inside the volume of the core. The de- 
nomination of one-neutron halo is somewhat misleading 
from that point of view. The influence (5it!haio of the halo 
on the nuclear extension is large, of about 0.4 fm out of a 
total root-mean-square (r.m.s.) radius of 2.90 fm. Note 
that such a value is very close to the value extracted ex- 
perimentally when going from ^^Be to ^°Be [53J. Last 
but not least, iVhaio and (5i?haio are found to be negligible 
for the two control cases considered, i.e. ^^C and the core 
density of ^^Be. This illustrates the ability of the method 
to discriminate between halo and non-halo systems. 



2. Medium-mass nuclei 

In Ref. [l|, the formation of neutron halos in Cr and 
Sn isotopes has been investigated through self-consistent 
EDF calculations performed using the SLy4 Skyrme func- 
tional complemented with a density-dependent mixed- 
type pairing functional (see Sec lIIIBTJ and Sec. lIIIB"2|) . 
The value of rg is determined using as an input the one- 
body density obtained from a symmetry-breaking HFB 
state, which is assumed to map out the nuclear internal 
density (see discussion in Ref. [l|). 

The (perturbative) excitation spectra of the last three 
bound odd Cr isotopes {Z = 24) shown in Ref. [l[ display 
optimal energy scales as far as the formation of a halo is 
concerned. The values of iVhaio and (5i?haio found for the 
last four bound even Cr isotopes, ^°Cr being the pre- 
dicted neutron drip-line nucleus for the parametrization 



used, are listed in Tab. |TT1 Beyond the A'' = 50 neutron 
shell closure, the steep increase of both Ny^hIo and (5i?haio 
is the signature of a halo formation. A decorrelated re- 
gion containing up to ^ 0.5 neutron appears at the limit 
of stability. Such a number of neutrons participating in 
the halo is small relatively to the size of the system but 
comparable in absolute value to the number found in light 
halo nuclei, as recalled in Sec. Ill Bl 

The contribution of the halo region to the nuclear ex- 
tension reaches about 0.13 fm in *°Cr. On the one hand, 
such a value is significant in comparison with the total 
neutron r.m.s. radius of those systems. On the other 
hand, it corresponds to about one-third of the value 
found for ^^Be. It is likely that the formation of halos is 
hindered as the mass increases because of the increased 
collectivity. We will come back to that in section |Vl 
Of course, one should not disregard the fact that single- 
reference EDF calculations miss certain long-range cor- 
relations which, once they are included, might slightly 
change the picture. 

One interesting feature discussed in Ref. [l| is that the 
kink of the r.m.s. radius at A'^ > 50, which had been in- 
terpreted as a halo signature in previous works, is partly 
due to a plain shell effect associated with the sudden 
drop of the two-neutron separation energy. As a matter 
of fact, the quantity (5i?haio allows the disentanglement 
of shell and halo effects in the increase of the neutron 
r.m.s. radius across N — 50. For isotopes further away 
from their drip-line, e.g. across TV = 82 in tin isotopes, 
the use of (5i?haio demonstrates that the kink of the neu- 
tron r.m.s. radius is entirely produced by the shell effect, 
whereas its further increase beyond A'^ = 82 is related to 
the growth of the neutron skin, not to the appearance of 
a halo [l|. Generally speaking the quantity SRhaio does 
not incorporate the contribution from the neutron skin 
and only characterizes the halo part of the density profile, 
e.g. spatially decorrelated neutrons. 





iVhalo 


(5-Rhaio [fm] 


^*Cr 


0.000 


0.000 


^«Cr 


0.057 


0.018 


^«Cr 


0.194 


0.055 


80 Cr 


0.472 


0.134 



Table II: Values of A^haio and SRhaio for the last bound even- 
even Cr isotopes, as predicted by the {SLy4-|-REG-M} func- 
tional. Results slightly differ from the ones shown in Paper I 
because of the different codes used to perform the EDF cal- 
culations. The difference remains very small, of the order of 
the third significative digit. 

A decomposition of the halo region in terms of canoni- 
cal states shows that it contains equal contributions from 
neutron 3si/2 and 2^5/2 orbitals. The possibility for 
£ >2 states to participate in the halo was not rejected 
from the outset. Eventually, the loosely bound 2^5/2 shell 
strongly contributes here because of its rather large occu- 
pation and intrinsic degeneracy. The latter observations 



point towards the formation of a collective halo in drip- 
line Cr isotopes, formed by an admixture of several over- 
lap functions, and call for a softening of the restrictions 
spelled out in light nuclei [0| . 

The energy scales necessary to the halo formation are 
not seen in the spectra of drip-line tin isotopes and 
the weak kink observed in the neutron r.m.s. radius 
at A^ = 82 is nothing but a shell effect, as mentioned 
above. In contrast with previous studies based on the 
Helm model [13, [H [H, [13, Ull , no collective halo is 
identified in tin isotopes, as the decorrelated region is 
found to have almost no influence on the matter exten- 
sion pj]. 



3. Atom-positron/ion-positronium complexes 

In atomic physics, valence electrons of neutral atoms 
can be located at large distances from the core. Because 
of the very long range of the Coulomb interaction, the 
penetration of the wave-function into the classically for- 
bidden part of the potential as the separation energy of 
the system becomes small cannot be interpreted as a halo 
formation ^6]. However, a positron can be attached to a 
neutral atom by the polarization potential, which can be 
parameterized as 



^pou(r) 
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(8) 
(9) 



where ap is the core polarization constant and f3 a cut- 



off distance. In this case, the 



decay of the po- 



tential at large distances does not ensure that parti- 
cles are able to tunnel through the potential barrier. It 
was found that several atorn-positron complexes can ex- 
ist [U, [H, M, M, M, M, M, [eaLand have been identified 
as having halo characteristics [12]. To quantify such an 
observation, values of A^haio and SRhaio are evaluated for 
such systems. 

The Hamiltonian of the atom-positron system with 
A^vai valence electrons reads with normalized units 
(me=e = l) ^^ 









-J 



-v^-Vdir(n)) + ypoi,i(n)) 



E R3^ - ^poi,2(rl, fo) , (10) 



\ri - ro 



where tq is the positron position vector, rij = ri — tj is 
the relative position of two valence electrons, whereas 
the direct Vdir and exchange T4x potentials between va- 
lence electrons and the core are computed exactly in the 
Hartree-Fock approximation. The two-body polarization 
potential is defined as 



Vpo\,2(fi,rj) 
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(11) 
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Figure 2: (Color Online) Valence electron, positron, core and 
total density for the e"'"-|-Be complex evaluated with the fixed- 
core stochastic variational method. All length units are nor- 
malized to the Bohr radius oq. The value of ro and the toler- 
ance margin are indicated by vertical lines. 
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Figure 3: (Color Online) Same as in Fig. [2] for the e"'"+Li 
system. 

When the system e++A is bound, its asymptotic be- 
havior can correspond to (i) a neutral core A plus a 
positron e"'", or (ii) a charged core A^ and a neutral 
positronium complex Ps, depending on the relative bind- 
ing energies of those configurations. Calculations are per- 
formed with the fixed-core stochastic variational Method 



(FCSVM) [iiliillSlIlllil, in a basis of explicitly cor- 
related gaussians for the individual wave-functions. The 
basis is taken large enough to correctly reproduce the 
asymptotic behavior of the e"'"-|-A or Ps-|-A+ systems. 
The results of such calculations for e+-|-Be, correspond- 
ing to the e^Be complex, are presented in Fig. [21 where 
the separation at large distances between a weakly bound 
positron and a core composed of the electrons is visible. 
As a result, a positron extended tail appears. In Fig. [3] 
are also displayed the results for the e^-l-Li system, which 
corresponds to a FsLi"*" complex. Indeed, one observes 
that the density tail is composed of almost identical e^ 
and e~ components. 

The results of the analysis, performed for several atom- 
positron complexes, are presented in Tab. IIIII The sep- 
aration energy e in the appropriate channel (e"'"-|-A or 
Ps+A+) is small compared to the ground state energy 
Egs of the complex. The situation regarding the energy 
scales at play is very favorable as far as the formation 
of halos is concerned. It is also possible to evaluate the 
composition of the halo region in terms of the proportion 
of electrons Pg- and of positrons Pe+ ■ 

The values of A^haio and i5i?haio demonstrate the ex- 
istence of halos in e^-l-Be, e"'"-|-Mg and e^-l-Cu which 
strongly affect the system extensions. For example, the 
spatially decorrelated part of the density accounts for 
about half of the total r.m.s. radius in e^-l-Be, although 
it contains only ^ 0.7 particle in average. In those cases, 
a positron halo is predicted as the halo region is al- 
most exclusively built from the positron wave-function 
(Pe+ ^^e-)- For e+-fLi and e++He, extremely large 
values for iVhaio and (5-Rhaio are extracted, demonstrat- 
ing that one is dealing in these cases with gigantic ion- 
positronium halos (Pe+ w Pe-)- 

Considering the values of iVhaio and (5Phaio, one real- 
izes that atom-positron and ion-positronium complexes 
display more extreme halo structures than nuclei. This 
is of course due to the nature of the interactions at play 
in such systems. 



C. Universality of the phenomenon 

As demonstrated, the method developed in Ref. [1| can 
be applied successfully to finite many-fermion systems of 
very different scales. This relates to the fact that the 
method relies on a model-independent analysis of the in- 
ternal one-body density P[i\[f). In all cases, a fraction 
of the constituents may extend far out from the core and 
influence strongly the size of the system. 

It happens that halo systems display scaling proper- 
ties which do not depend on their dimension and con- 
stituency. This can be characterized by the extension 
of the halo wave- function {r\) as a function of the sep- 
aration energy E, and those quantities can be made di- 
mensionless using as a scale the classical turning point 
Rq for the interaction potential of interest and the re- 
duced mass of the systems /x 0, [l^, [3J, [33| ■ The generic 



Atom 


Asympt. 


N,- 


A^halo 


-Rr.m.s. [do] 


SRhalo [fflo] 


n+ [%] 


P.- [%] 


Egs [at. units] 


e [at. units] 


Ref. 


Be 


e++Be 


4 


0.624 


5.661 


3.194 


98.1 


01.9 


-1.0151 


0.0032 


[12] 


Mg 


e++Mg 


12 


0.669 


2.298 


0.826 


80.3 


19.7 


-0.8477 


0.0156 


[52] 


Cu 


e++Cu 


29 


0.754 


1.777 


0.975 


88.6 


11.4 


-0.2891 


0.0051 


m 


He 


Ps+He+ 


2 


1.982 


15.472 


14.568 


50.3 


49.7 


-2.2506 


0.0006 


[66] 


Li 


Ps+Li+ 


3 


1.972 


7.781 


7.088 


50.8 


49.2 


-7.5324 


0.0024 


[67] 



Table III: Results of the halo analysis for various atom-positron systems evaluated with the fixed-core stochastic variational 
method. The columns provide: neutral atom symbol, asymptotic form, total number of electrons A'^^.- , halo factors Ahaio and 
SRhaio, total matter r.m.s. radius -Rr.m.s., relative proportions P^+ and P^- of electrons and positrons in the halo region, ground 
state energy -Bgs of the e"''-|-A complex, and binding energy e with respect to the corresponding dissociation threshold. All 
length units are normalized to the Bohr radius ao, and energy values are in atomic units (1 at. unit= 27.21162 eV). 



asymptotic scaling laws of the two-body system are ex- 
tracted using a finite spherical well, and depend on the 
angular momentum of the weakly-bound overlap func- 
tion [SJ], as seen in Fig. S) Results for Hght nuclei obey 
rather well such universal scaling laws. For few-body sys- 
tems, it is commonly admitted that halos appear when 
{r\) / R^ > 2, which corresponds to a probability greater 
than 50 % for the weakly-bound nucleon to be in the for- 
bidden region [y]. 

The results obtained for medium-mass nuclei can also 
be displayed in Fig. [4] and compared with the generic 
scaling laws. However, the dimensionless quantities have 
to be redefined. For medium-mass systems the halo r.m.s. 
radius is evaluated, by analogy with Ref. [35], through 
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(12) 



where {r^^t) is the total neutron r.m.s. radius, while the 
core r.m.s. radius is approximated by -Rr.m.s., inner (see 
Eq. ([3|)). The reduced mass is taken as the effective 
isoscalar nucleon mass m*, while the classical turning 
point Rq of the central part Uq (r) of the one-body poten- 
tial is evaluated, by analogy with the finite- well potential, 
as fH 



drr^Uq{r) 
dr r Uq (r) 



Rl 



(13) 



where q — n as we are interested in neutron halos. 

In fig. m the last bound Cr isotopes are located in- 
between the ^ = and i = 2 scaling curves, with ^^Cr 
being closer to the i = 2 curve than *°Cr. This is con- 
sistent with the admixture of orbitals that builds the 
corresponding halos, as discussed in Sec. Ill B 21 Ij. The 
neutron density of most medium-mass halo nuclei does 
not extend as much as those of few-body systems such 
as "Be. Still, '^^Cr and ^°Cr display few-body-hke halo 
properties and the ratio {r\) / R^ does exceed 2 for ^"Cr. 
On the contrary, the extension of neutron-rich tin iso- 
topes is not significant enough in regard with their sepa- 
ration energy to be characterized as halo systems. This 
is consistent with the findings of Paper L 



It would be of interest to place atom-positron com- 
plexes in Fig.m However, the many-body method used to 
compute their properties does not allow an easy extrac- 
tion of the corresponding classical turning point. As a 
result, the corresponding results do not appear in Fig. 3] 
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Figure 4: (Color Online) Universality of halos features. 
The dimensionless extension of the halo is plotted against 
the dimensionless separation energy for EDF calculations of 
medium-mass nuclei (filled symbols) and experimental results 
for few-body systems (open symbols). Generic scaling rela- 
tionships obtained from a finite spherical well are given in 
solid lines, while the gray-shaded area corresponds to the halo 
are, as defined by the commonly-used criterion for few-body 
halo nuclei {rl)/R^ > 2. 



III. TECHNICAL ASPECTS 

The objectives of the remaining part of the present 
study are (i) to predict which spherical medium-mass 
nuclei might display halo features, (ii) to check the de- 
pendence of the results on several ingredients of the nu- 
merical implementation, (iii) to probe the sensitivity of 
the predictions to the characteristics of the many-body 
treatment and (iv) to study the specific impact of pairing 
correlations on halo nuclei. 



In the present section, and because halos are extreme 
systems whose asymptotic must be properly accounted 
for, we study the dependence of our predictions on certain 
technical features of the many-body calculation. To do 
so, ingredients of single-reference EDF calculations are 
briefly recalled at first. 



B. Energy density functional 

1. Particle-hole channel 

The Skyrme part of the EDF takes, in the case of spher- 
ical nuclei, the standard form [321 



A. HFB equations 



£PP = —TO + V CP.pl + C^PpT Apt + C^pTTT 

}Tn ^ — ^ 



T=oa 



X C^tJt + C^ PT^ ■ Jt , (15) 



Spherical symmetry is assumed throughout the rest of 
the paper and spin/isospin indices are sometimes omitted 
for simplicity. Calculations are performed using a code 
that takes advantage of the so called "two-basis method" 
to solve the HFB equations [68| . Thanks to the spherical 
symmetry, the HFB equation are solved for each (/,j) 
block separately in the basis {4>nij} that diagonalizes the 
single-particle field h 
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where f/^^ and Vll 



(14) 
are the expansion coefficients of the 
upper and lower parts of the HFB spinor ^nij on the 
basis {\(t)klj)}^ whereas h denotes the pairing field. The 
two-basis method authorizes to perform calculations in 
very large boxes. The nucleus is put in a spherical box 
such that wave functions are computed up to a radial dis- 
tance i?box, with vanishing boundary conditions (Dirich- 
let) imposed. The value of i?box has to be chosen large 
enough as to ensure convergence of the calculations (see 
Sec. lHIDJ) . The differential equation to find the {(t^nij} is 
solved on a discrete mesh of step size h = 0.25 fm using 
the Numerov algorithm [6^, UM ■ 



In Eq. (|14p the upper bound of the sum over the index 
k is not specified. In an actual calculation, the sum is 
truncated by keeping states {4>nij} up to a certain max- 
imum energy i?max- Its actual value ranges from several 
MeV up to hundreds of MeVs depending on the method 
used to tackle the ultra-violet divergence of the local 
pairing functional (see section IIIIB ip . In addition to 
the energy cut, a truncation j'^^ on the number of par- 
tial waves kept in the basis {ipnij} is implemented. In 
principle, all wave functions below £',nax should be kept, 
but this makes the computation time rather long whereas 
the wave functions with very high angular momenta will 
not contribute to the nuclear density. Nonetheless, such 
a truncation must not be too drastic, in particular for 
loosely bound nuclei. Checks of convergence with respect 
to -Rbox, -Emax and j'^ut of observables of interest in halo 
nuclei are discussed in Sec. IIIIDI 



where p{f)^ t(j^) and J(r) denote the normal, kinetic 
and spin-orbit densities, respectively. The parameter x is 
equal to zero or one, depending on whether tensor terms 
are included or not in the functional, while the index T 
labels isoscalar (T — 0) and isovector (T — 1) densities. 
For protons, the particle-hole part of the EDF is com- 
plemented with a Coulomb term whose exchange part is 
treated within the Slater approximation |32l | . 

To study the effect of specific features of the particle- 
hole functional on the formation of halos, a set of Skyrme 
functionals characterized by different properties is used 
in the present study: (i) SLy4 stands as a reference 
point, (ii) SHI displays a different density dependence 
which leads to a too high infinite matter incompress- 
ibility i^oo, (iii) T6 has an isoscalar effective nucleon 
mass {m*/m)s = 1, providing a denser single-particle 
spectrum, (iv) SKa has a low isoscalar effective mass 
and a different density dependence (density-dependent 
term with an exponent of 1/3 instead of 1/6), (v) the 
functional "to*1" has been specifically adjusted for the 
present work with the same procedure as for SLy4 but 
with the constraint (m*/m)s = 1, (vi) the parameteriza- 
tions "Psat " have also been adjusted specifically, with 
different nuclear matter saturation densities Psat, (vii) 
T21 to T26 incorporate tensor terms that differ by their 
neutron-neutron couplings [3^. The parameterizations 

1 /2 /3 

"m*l" and "yOaat " have been adjusted using the proce- 
dure of Ref. [39| which amounts to reproducing (i) the 
binding energies and charge radii of ^°''*^Ca, ^^Ni, ^°Zr, 
^^^Sn and ^°*Pb, (n) the binding energy of ^°°Sii^ and 
(iii) the equation of state of pure neutron matter [7l| and 
other standard properties of symmetric nuclear matter as 
well as the Thomas-Reiche-Kuhn enhancement factor of 
the isovector giant dipole resonance. 

Infinite matter properties of all used parameterizations 
are summarized in Tab. lIVl The isovector effective mass 
(related to the Thomas-Rciche-Kuhn enhancement fac- 
tor Ky) is significantly different for these parameteriza- 
tions but its effect on static properties of nuclei is rather 
smaU U^. 



2. Particle-particle channel 
The local pairing functional used 
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K^ 


{m*/m)s 


K« 


E/A 


Ref. 


SLy4 


0.160 


229.9 


0.70 


0.25 


-15.97 


[73,74J 


SIII 


0.145 


355.4 


0.76 


0.53 


-15.85 


[75j 


m*l 


0.162 


230.0 


1.00 


0.25 


-16.07 


[76] 


pLt 


0.145 


230.0 


0.70 


0.25 


-15.69 


[76] 


pL 


0.160 


230.0 


0.70 


0.25 


-15.99 


[76] 


pL 


0.175 


230.0 


0.70 


0.25 


-16.22 


[76] 


T6 


0.161 


235.6 


1.00 


0.00 


-15.93 


[77] 


SKa 


0.155 


263.1 


0.61 


0.94 


-15.99 


[78] 


T21-T26 


0.161 


230.0 


0.70 


0.25 


-16.00 


[39] 



Table IV: Infinite matter properties of Skyrme functionals 
used in the present study: saturation density psat [fm~^], bulk 
compressibility Koo [MeV], isoscalar effective mass {m* /m)s, 
Thomas-Reiche-Kuhn enhancement factor Ki, and energy per 
particle at saturation E/A [MeV]. 
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derives from a Density-Dependent Delta Interaction 
(DDDI) with the same strength VO for neutron-neutron 
and proton-proton pairing. In addition, two parameters 
ry and a control the spatial dependence of the effective 
coupling constant. With psat designating the saturation 
density of infinite nuclear matter, a zero value of rj corre- 
sponds to a pairing strength that is uniform over the nu- 
clear volume ( "volume pairing" ) while t] — 1 corresponds 
to pairing strength which is stronger in the vicinity of 
the nuclear surface ("surface pairing"). A value ry = 1/2 
corresponds to an intermediate situation ("mixed-type 
pairing"). The parameter a is usually set to one. Values 
a < 1 correspond to stronger pairing correlations at low 
density. In the present work, we are interested in varying 
such empirical parameters over a large interval of values 
to quantify how much the characteristics of the pairing 
functional impact halo systems. Note finally that, all 
along this work, the strength Vq is chosen so that the neu- 
tron spectral gap (A)?) [79] equals 1.250 MeV for ^^OSn. 
Such a value of Vq provides reasonable gaps in Ca, Sn 
and Pb regions. 

To compensate for the ultra-violet divergence of the 
pairing density generated by the local pairing functional, 
a common procedure consists of re gula rizing all integrals 
at play through the use of a cutoff [80|, e.g. on quasipar- 
ticle energies Ei < Ecut- Pairing functionals using such 
a regularization scheme are referred to as "REG" in the 
following. In particular REG-S, REG-M and REG-V de- 
note regularized surface-, mixed- and volume- type pair- 
ing functionals, respectively. If the parameter 77 differs 
from 0, 1/2 and 1 or if a differs from 1 the functional is 
generically noted as REG-X. 

Using such a regularization method, the pairing 
strength is adjusted for each cutoff E^ut, the latter be- 
ing eventually taken large enough for observables to be 



insensitive to itsprecise value. A widely used value is 
^cut = 60 MeV [7i,[83,[8l|. As the density dependence of 
the pairing functional is made more surface-peaked, the 
pairing strength increases, for a fixed value of Scut, as is 
exemplified in Fig. [5l This is a consequence of the fitting 
procedure that uses a single nucleus to adjust the overall 
strength. Indeed, if the pairing strength is peaked at the 
nuclear surface, individual gaps decrease, especially for 
well-bound orbitals residing in the nuclear interior. To 
compensate for this effect and maintain the same value 
of (A") in one given nucleus, the overall pairing strength 
has to be increased. 




Figure 5: (Color Online) Pairing strength Vb as a function of 
77 and a for regularized REG-X and in-medium renormalized 
REN-X pairing functionals. 



As an alternative to the sharp cut-off regularization, 
one can identify and throw away the diverging part of the 
pairing density through the use of an auxiliary quantity 
that diverges in the same way but from which the diver- 
gence can be subtracted analytically We follow here the 
procedure introduced in Refs. [82, [83, Isl, [S^, ISal • Using 
such a renormalization procedure, results become inde- 
pendent of Ecut reaches about a few tens of MeVs [85| . 
In the following, calculations of finite nuclei are per- 
formed using a conservative value of Ecut = 200 MeV 
(see Sec. lIIIDSt . Pairing functionals combined with the 
renormalization scheme are referred to as "REN" in the 
following. In particular REN-S, REN-M and REN-V 
denote renormalized surface-, mixed- and volume- type 
functionals, respectively. If the parameter 77 differs from 
0,1/2 and 1 or if a differs from 1 the functional is gener- 
ically noted as REN-X. 



C. Regularization versus renormalization scheme 

It was observed in Ref. [83] that combining the sharp 
cut-off regularization method with an extreme surface 
pairing functional (a < 1) might lead to unrealistic mat- 
ter and pairing densities. As our goal is to study the 



influence of the pairing functional attributes on halo 
properties, it is of importance to characterize such a 
feature in more detail. As a matter of fact, it can 
be shown that extreme surface REG-X functionals lead 
to a spurious Bardeen-Cooper-Schrieffer (BCS) Bosc- 
Einstein-Condensation (BEG) phase transition of infi- 
nite matter. On the contrary such a nonphysical tran- 
sition, associated with the use of a pairing functional 
that (wrongly) predicts a bound di-neutron system in 
the zero-density limit [8g |. does not occur with extreme 
surface REN-X functionals [8a |. 

In the present section, we briefly explain that such 
a spurious feature is in one-to-one correspondence 
with producing unphysical matter and pairing densities 
through HFB calculations of finite nuclei. To do so, prop- 
erties of the last bound Gr isotopes are calculated using 
different pairing functionals (i) (77 < l,a ~ 1), which 
covers from volume to standard surface pairing function- 
als, and (ii) (ry = 1, a < 1), which correspond to extreme 
surface pairing functionals. 

When the regularization scheme is used, one first ob- 
serves that the position of the drip-line seems to change 
drastically for extreme surface pairing functionals. For 
instance, ®'*Gr is predicted to be bound against neutron 
emission (A' < 0) [SOJ for a = 0.5. However, the situa- 
tion is more subtle than it looks at first. As exemplified 
by Fig. [HI for ^"Gr, nuclei are in fact no longer bound as 
a gas of low-density superfluid neutron matter develops 
as a decreases. The normal and pairing (neutron) densi- 
ties grow at long distances and become uniform beyond 
a radius r « 9 fm for a < 0.5 which is precisely the 
critical value for which infinite nuclear matter undergoes 
an unphysical BGS-BEG transition f89l|. It becomes in 
fact energetically favorable for the nucleus to drip bound 
di-neutrons and create a superfiuid gas. 

Two remarks are at play at this point. First, it is im- 
portant to use enough partial waves in the single-particle 
basis to describe the gas properly. Otherwise, one may 
not resolve it and conclude that the nucleus displays un- 
usually extended normal and pairing densities [87| . Here, 
calculations are performed including all partial waves up 
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151/2. This is of course an ex- 



treme situation and less partial waves, but not too few, 
have to be used to describe genuine halos as is explained 
in the next section. Second, the appearance of a uni- 
form gas of bound di-neutrons is driven by the neutron- 
neutron pairing interaction and not by the proximity of 
the continuum. As a matter of fact, the fictitious nu- 
cleus whose densities are displayed on Fig. [HI and which 
is held by the box, is still bound against singie-nucleon 
emission (A^ < 0). If we were to increase i?box, the sys- 
tem would turn into a gas of di-neutrons with binding 
energy 2A'' |88j . 

Finally, one concludes that using a functional which 
wrongly predicts the existence of a bound di-nucleon 
state in the ^^o channel in the vacuum translates into 
the creation of a spurious low-density di-neutron gas sur- 
rounding finite nuclei. Such an observation is of impor- 
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Figure 6; (Color Online) Neutron normal p"(r) and pairing 
p"(r) densities plotted in logarithmic scale for *"Cr and sev- 
eral regularized REG-X pairing functionals. The calculations 
are done in a box of 40 fm and using all partial waves up to 

.?cut ^ .7cut ^ <^max ^ iOl/2. 




Figure 7: (Color Online) A^haio and 5_Rhaio parameters for Cr 
isotopes and different regularized REG-X pairing functionals. 



tance regarding the analysis of halo systems as it signals 
that the use of strongly surface-peaked pairing function- 
als, combined with the regularized scheme, might lead 
to un-physical predictions. The danger resides in par- 
ticular in the use of pairing functionals which are not 
obviously un-physical, i.e. for which the di-neutron gas 
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is not yet fully developed*^^). In such cases, the calcula- 
tion will lead to wrongly predict the existence of gigantic 
halos, e.g. see in Fig. [7] the unreasonable values of A^haio 
and SRhaio predicted for the last bound Cr isotopes by 
REG-X functionals with a < 1. 

From that point of view, the renornialization scheme 
is safer as it prevents the problem from happening. In- 
deed, normal and pair neutron densities remain localized 
and evolve very little as a < 1 decreases [80|. On the 
other hand, it is crucial to point out that no problem 
occurs with the regularization scheme cither as long as 
the low-density part of the pairing functional is physi- 
cally constrained, e.g. (77, a) are adjusted as to reproduce 
pairing gaps calculated in infinite nuclear matter through 
ab-initio calculations [90, l9l| ■ Eventually, the standard 
fitting strategy used here is the real source of the poten- 
tial problem rather than the regularization method itself. 
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D. Convergence checks 

Several basis truncations are introduced under the 
form of (i) a box of finite radius i?box, (h) an angular- 
momentum cutoff j^^^ for each isospin and (iii) a contin- 
uum energy cutoff i^max, to accelerate the convergence 
of the calculations. Such truncations are physically mo- 
tivated as (i) the nucleus is localized in space (ii) high- 
lying unoccupied states are not expected to contribute 
to ground state properties. However, the values of the 
truncation parameters have to be carefully chosen not to 
exclude meaningful physics. As a result, the convergence 
of observables of interest has to be checked, which we do 
now. All tests are performed for the chromium isotopes 
using the {SLy4-|-REN-M} functional. 



1. Box radius 

The evolution of the halo factors, which are the 
most critical observables related to halo properties 
are represented in Figs. [8] for *°Cr as a function of 
-Rbox, for angular-momentum truncations j^^^ = 65/2 
and j^yt = 61/2, and an energy cutoff ii^max = 200 MeV. 
The pairing strength is not refitted for each bin, since 
the overall effect is found to be negligible. For small 
box sizes, iVhaio and (5i?haio are not fully converged and 
increase with i?box, SRhedo being the most sensitive quan- 
tity. A box radius i?box — 40 fm, for which convergence 
is achieved for all observables of interest, is used in the 
following. 



Figure 8: (Color Online) Evolution of A^haio and iS-Rhaio as a 
function of the box radius for *'^Cr. 



2. Angular-momentum truncation 

The choice of j"^^ and j^^^. is critical because the 
partial-wave truncation impacts the way the continuum 
is represented in the calculations. Indeed, high-^ orbitals 
can contribute to nuclear properties such as the pairing 
gaps. To quantify such an effect, we calculate, in the 
canonical basis, the probability distribution of particle 
P'^aO) and pair -P^^(j) occupations as a function of the 
single-particle angular-momentum j 
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^ Most probably because of inappropriate numerical parameters. 



The neutron distribution P^\ (j) is shown in Fig. [9] for 
the halo nucleus ^''Cr and different angular momentum 
truncations j"^j = j^^^ = Jmax- It is found that the v'^'^ 
strength is mostly distributed over states with j < 9/2. 
As a result, such an occupation distribution is converged, 
at least to first approximation, for j"^j — j^^^ — Jmax = 
15/2. The corresponding P^y{j) distribution (Fig. [TU)) 
extends much further towards high j values. This could 
be expected as w' is maximum for deeply bound canoni- 
cal states whereas m' i;^ is peaked around the Fermi level 
and decay slower as the canonical energy increases above 
A'. Correspondingly, the local pair density p''(r) extends 
further out in space than the normal density p'^(r) [80]. 
As one goes to drip-line (halo) nuclei in particular, p'^(r) 
extends very far out and requires many partial waves 
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to be well represented. This was already highlighted in 
Ref. [931 . It is clear from Fig. \lU\ that a minimal cutoff 
of Jmax = 31/2 is needed to achieve a reliable enough 
description of the pair distribution. One sees in particu- 
lar that for Jmax — 15/2 (i) the missing u" w" strength at 
high j is wrongly redistributed over j — 1/2 and j — 5/2, 
making those states more paired and thus more localized, 
whereas (ii) some of the f "^ strength of the j = 1/2 states 
is transferred to j = 5/2 states. Considering that 3si/2 
and 2(i5/2 states are precisely those building up the halo 
in ^"^Cr, the latter two effects associated with using a two 
small angular momentum cutoff inhibits artificially the 
formation of the halo. The results are qualitatively iden- 
tical for the regularization scheme: in both cases, rather 
high values of Jmax are needed to properly describe the 
continuum. 



0.4 



0.3 



C ^0.2 



0.1 



0.0 



Jn.ax=312 



80 



Cr 



13 19 25 31 37 43 



2j 



Figure 9: (Color Online) Distribution of neutron canonical 
states occupations «" ^ as fractions of the total strength for 
*°Cr computed with {SLy44-REN-M} functionals and differ- 



ent angular momentum truncations j"^, 
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Figure 10: (Color Online) Same as Fig. [9] for the distribution 

of M"t)". 

Indeed, this is what is observed in Fig. II II where A^haio 



and SRhaio are given for ®°Cr as a function of Jmax- 
quantities reach a plateau around Jmax = 31/2. 
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Figure 11: (Color Online) Evolution of A'l^aio and (5i?haio in 
**°Cr and {SLy4-fREN-M} functionals for different angular 
momentum truncations jSit = Jcut ~ Jmax • 



The previous analysis shows that HFB calculations 
with too small values of j"^^ and j^^j cannot be trusted 
at the limits of stability if one is interested in detailed in- 
formation about potential halos. Of course, considering 
the ultimate experimental accuracy achievable for mat- 
ter r.m.s., one should not be too extreme as far as the 
required convergence is concerned. From a theoretical 
perspective, and considering the theoretical error bars 
on the determination of Miaio and 5 JJhaio , it is necessary 
to include partial waves up to j^^^ 



Jcut ~ Jmax — ol/z. 



3. Energy cutoff 

The value of the energy cutoff -Emax in the quasiparticle 
continuum is an important parameter of the calculation. 

For regularized pairing functionals, the values of Jmax 
and i^max must be taken large enough that, including 
a renormalization of the coupling strength through the 
re-fitting of data, the observables of interest do not de- 
pend on their precise values. It was found that smaller 
basis truncations could be used for the REG case than 
for the REN case, as convergence is reached faster, as 
exemplified in Figs. [12] and [131 For the REG case, con- 
vergence for the ground state energy as well as for the 
neutron pairing gap is almost achieved for Jmax = 21/2 
and E^max = 200 MeV. 

For renormalized functionals, the situation is more sub- 
tle. First, Emax must be large enough for the result to 
be independent of its value [83|. However, it must be re- 
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membered that the field theory renormahzation scheme 
subtracts a diverging part on the basis that aU partial 
waves below a certain energy cutoff have been included. 
Thus, for a given (high enough) -Emax, the angular mo- 
mentum truncation must be large enough to prevent the 
counter term from removing contributions of states that 
were not considered in the first place. This is illustrated 
in Figs. [HI and 1151 which display the binding energy and 
neutron spectral gap of ^"Cr as a function of the angular 
momentum cutoff Jmaxj for fixed values of i?max- Note 
that all ii^niax values considered are large enough to obtain 
converged observables. However, one sees that increasing 
the energy cutoff necessitates a larger number of partial 
waves to reach the converged values for both the energy 
and the gap. Consequently, it can be counter-productive 
to use a safe but unnecessarily large energy cut -Emax 
as it results in the necessity to also increase Jmax- On 
the other hand, the proper description of certain physi- 
cal phenomena such as halos intrinsically requires a large 
number of partial waves as discussed in the previous sec- 
tion. In such a case, one first fixes Jmax and makes sure 
to use a coherent energy cutoff. In the present work, 
we use -Emax = 200 MeV and Jmax = 65/2. This corre- 



sponds to an conservative choice as E„ 
Jmax = 41/2 would be sufficient. 
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Figure 13: (Color Online) Same as in Fig. ll2l for the neutron 
spectral gap. 
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Figure 14: (Color Online) Same as in Fig. [12] for the 
{SLy4-FREN-M} functional. 
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Figure 12: (Color Online) Binding energy as a function of the 
angular momentum cutoff Jmax for fixed values of Emax as 
obtained from {SLy4+REG-M} functional for *°Cr. 



IV. IMPACT OF PAIRING CORRELATIONS 

A. Pairing anti-halo effect 

In the presence of pairing correlations, the asymptotic 
of the one-body neutron density takes a different form 
from the one it has in the EDF treatment based on an 
auxiliary Slater determinant [l|. Indeed, the decay con- 
■Cg = v^— 2meg/?i2, with |e 



slant reads as 



<r9| — 
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Eq = Tain^lE^] being the lowest quasiparticle excitation 



0.50 



'Emax=70 MeV 
'En,a,^=200 MeV 
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'E„a,= 1000MeV 
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Figure 15: (Color Online) Same as in Fig. [13] for the 
{SLy4-hREN-M} functional. 
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energy. Considering a canonical state 0o lying ^t the 
Fermi level near the drip- line (bq w A' « 0), one finds 
that Eq w Aq > 0. Therefore, in first approximation, 
paired densities decrease faster than unpaired ones and 
pairing correlations induce an anti-halo effect by localiz- 
ing the one-body density OEHIii- 

To evaluate the quantitative impact of this effect, drip- 
line Cr isotopes have been calculated with and without 
explicit treatment of pairing correlations. In the latter 
case, a filling approximation [9J] is used for incomplete 
spherical shells. In both cases, the SLy4 Skyrme func- 
tional is used. When including pairing, a mixed-type 
pairing is added. A comparison between neutrons single- 
particle levels is represented in Fig. [16] for the last bound 
nuclei, ^^Cr being predicted to be bound when pairing 
correlations are excluded from the treatment. This is 
interesting in itself as it shows that pairing correlations 
can change the position of the drip-line and modify in 
this way the number of halo candidates over the nuclear 
chart. There is only little difference between the canon- 
ical energies e" in the two cases^^^. However, the values 
of the halo criteria iVhaio and 5i?haio are significantly dif- 
ferent, i.e. the neutron halo is significantly quenched in 
®'^Cr when pairing is omitted whereas the situation is 
reversed in the lighter isotopes, as seen in Fig. [181 
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Figure 16: (Color Online) Neutron single-particle energies e" 
along the Cr isotopic chain. Left: eigenvalues of h" obtained 
when omitting pairing correlations in the calculation. The 
non-resonant continuum spectrum has been omitted while 
resonant states with positive e" are drawn as discrete states 
for convenience. Right: canonical energies obtained including 
pairing correlations through a mixed-type REG-M functional. 
The Skyrme SLy4 functional is used in both cases. The con- 
ventions used for single-particle states labeling are identical 
to Ref. 1| and are recalled in Fig. [TTl The Fermi level (ffl) is 
set to the last filled orbital for the EDF treatment based on 
an auxiliary Slater determinant. 
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Figure 17: (Color Online) Conventions used in all the fig- 
ures for the labeling of individual states and of the chemical 
potential. 

Such results underline that pairing correlations affect 
halos in two opposite ways. Pairing (i) inhibits the forma- 
tion of halos through the anti-halo effect (ii) enhance the 
formation of halos by scattering nucleons to less bound 
states with smaller decay constants. For example, the 
anti-halo effect dominates in ^^Cr and ^*Cr whereas the 
promotion of neutrons into the weakly bound 3si /2 makes 
the halo to be more pronounced in ^"^Cr when pairing is 
included. 

When pairing is omitted, the number of nucleons in 
the 2dt^/2 orbital increases linearly as one goes from ^^Cr 



to 



80 



Cr. At the same time, the 2^5/2 shell becomes more 
bound as a result of self-consistency effects. This explains 
the negative curvature of (5i?haio as one goes from '^^Cr 
to ^"Cr, while A'haio is indeed almost linear. In ^^Cr, 
two effects contribute to the very pronounced halo that 
is predicted in the absence of pairing correlations (i) the 
3si/2 gets fully occupied whereas (ii) the pairing anti- 
halo effect is inoperative. The results discussed above 
are further illustrated in Fig. [TO] where the contributions 
of different single-particle states to the halo are shown. 

To specifically extract the pairing anti-halo effect, we 
now perform a toy model calculation of a fictitious ^°Cr* 
nucleus. Filling the single-particle wave- functions ob- 
tained from the calculation of *^Cr without pairing cor- 
relations with the canonical occupations obtained from 
the calculation of ^"Cr with pairing, we extract A'haio,! 
from each orbital*^ ■^•'. Such a procedure allows one to iso- 
late, in semi-quantitative manner, the net change of iVhaio 
due to the difference in the asymptotic of paired orbitals 
keeping their occupation fixed. Doing so, it is found that 
the contribution of the 3si/2 state to the halo is smaller 



^ The canonical basis is identical to the oigenbasis of h is the zero- 
pairing limit. 



^ Single-particle states extracted from the calculation of ^''Cr with- 
out pairing cannot be used because the essential 3si/2 orbital 
belongs to the continuum and has plane wave asymptotic in this 
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Figure 18: (Color Online) Halo criteria A'haio and SR^aio 
for Cr isotopes, obtained through spherical HF (solid lines) 
and HFB (dashed lines) calculations with the Skyrme SLy4 
functional. For HFB calculations, mixed-type REG-M pairing 
is used. 
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Figure 19: (Color Online) Contributions of individual orbital 
to the halo: EDF calculation of the least bound Cr isotopes 
with (bottom) and without (top) pairing. The 2di,/2 orbital 
corresponds to vertically-dashed areas, and 3si/2 orbital to 
filled ones. The other contributions are found to be negligible. 



in 8"Cr than in ^"Cr* by about 26%. For the 2^5/2 or- 
bital the suppression of A^haio.i is about 15% (the value 
of To is also different for *°Cr and ^"Cr*, but this does 
not affect the results significantly). As a result, one sees 
that if the anti-halo effect were ineffective, the scattering 
of particles into higher-lying orbitals would bring A^haio 
from 0.300 to 0.582. Instead, iVhaio is only increased to 
0.473 in the full fledged calculation of ^^Cr with pairing 
correlations, i.e. the anti-halo effect reduces the potential 
increase by 40%. 

To conclude, there is no simple answer as to whether 
pairing correlations enhance or hinder the formation of 
halos. The net result depends on structure details of the 
particular nucleus of interest JQ^, l9y| . 



B. Decorrelation from the pairing field 

An additional effect might come into play as far as the 
role of pairing in the formation of halos is concerned. 
Very weakly bound i — orbitals are expected to de- 
couple from the pairing field as the neutron separation 
energy goes to zero (42, |43, |44|- ^^ ^ result, such an 
orbital would not be subject to the anti-halo effect and 
may develop a very long tail. 

The signature of this effect can be seen in the single- 
particle occupation profile. The canonical occupations of 
all neutron single-particle states in all drip-line Cr iso- 
topes are gathered in Fig. [^D] and plotted as a function 
of (e" — A"), where e" is the neutron canonical energy 
and A" the neutron Fermi level. Those occupations are 



compared to the BCS formula 



,,n2 



v/(er - A")2 + {A^J^ 



(19) 



calculated using the maximum/minimum neutron spec- 
tral gap (A") found among all drip-line Cr iso- 
topes. The calculations have been performed using the 
{SLY4-fREG-M} functional. 

The s-wave occupation probability follows closely the 
BCS-type profile calculated using the minimal spectral 
gap. The high-^ orbitals follow well the BCS-type profile 
computed with the large spectral gap. This corroborates 
the trend discussed in Refs. [42, H, |4J] and underlines 
that £ ^ orbitals with e" ~ A" ~ are less paired 
than high-^ ones. This is also confirmed by looking at 
individual gaps given in Fig. 1211 The 3si/2 state displays 
a smaller gap than other orbitals as it approaches the 
Fermi level (from above), the latter reaching the contin- 
uum threshold. 

On the other hand, the canonical gap of the 3si/2 state 
remains significant as it crosses the Fermi level and the 
anti-halo effect is still in effect for that orbital. In fact, 
the critical decoupling of s orbitals from the pairing field, 
discussed in Refs. [i^, |43, [4J] through schematic HFB 
calculations, becomes fully effective only when the Fermi 
level and the s state energy are both of the order of a 
few keVs. Such extreme situations (i) are not reached in 
realistic systems [97[ (ii) would require an accuracy on 
the predicted value of the separation energy which is far 
beyond present capacities of EDF calculations. 
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**"& (pairing) 


iVhalo 


0.473 




e? [MeV] 


V?2 


NhBlo,i 


3Sl/2 
2d5/2 

159/2 

1/5/2 

2pi/2 

Other 


>0.0 
-0.173 
-0.665 
-4.056 
-8.673 
-8.920 
< -10.0 


0.451 
0.828 
0.993 
0.999 
0.999 


~ 1.3.10^2 
0.225 
0.222 
0.008 
0.001 
0.002 
~ 2.8 X 10"* 




**^Cr (no pairing) 


A^halo 


1.007 




el [MeV] 


VT2 


-/Vhalo.i 


3Sl/2 
2d5/2 
I59/2 
1/5/2 
2pi/2 

Other 


>0.0 
-0.312 
-0.858 
-4.200 
-8.783 
-9.078 
< -10.0 


1.000 
1.000 
1.000 
1.000 
1.000 


0.000 
0.675 
0.314 
0.010 
0.001 
0.002 
~2.7 X 10-^ 




^^Cr* (pairing) 


iVhalo 


0.582 




e? [MeV] 


<2 


A'halo.i 


3Sl/2 
2d5/2 
I59/2 
l,f5/2 
2pi/2 

Other 


>0.0 
-0.312 
-0.858 
-4.200 
-8.784 
-9.078 
< -10.0 


0.451 
0.828 
0.993 
0.999 
0.999 


?????? 

0.305 
0.260 
0.010 
0.002 
0.002 
~ 2.7 X 10"^ 



Table V: A^haio, individual contributions to the halo A^haio,! 
(see Paper I for the definition of A^haio,!) 1 single-particle canon- 
ical energies e" and occupations v" for (i) Cr with pairing 
(ii) ®^Cr without pairing, and (iii) ^''Cr* with pairing but us- 
ing the single particle spectrum calculated in the absence of 
pairing. 



C. Importance of Ioav densities 

The values of rj and a entering the pairing functional 
strongly affect the spatial localization of the pairing field, 
and thus the gaps of weakly-bound orbitals lying at the 
nuclear surface. 

In previous studies 87], it was found that the size of 
a neutron halo could change by one order of magnitude 
when the pairing functional evolves from a volume to a 
extreme-surface type. However (i) the evaluation of the 
halo size was performed through the Helm model, the 
limitations of which have been pointed out in Ref. [l|, 
and (ii) the standard regularization scheme was used with 
extreme-surface pairing functionals which, as discussed 
in Sec. IHICl leads to unsafe predictions of halo proper- 
ties. The renormalization scheme is used in the present 
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Figure 20: (Color Online) Neutron canonical occupations ii" 
as a function of (e" — A") for all drip-line Cr isotopes. The 
calculations are done with the {SLY4+REG-M} functional. 
The profiles computed from the BCS-type formula using 
the minimum (dashed line) and maximum (dashed-dotted 
line) spectral neutron gaps (A") among all those isotopes 
are shown for reference. Conventions for labeling individual 
states are given in Fig. 1171 
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Figure 21: (Color Online) Individual neutron canonical gaps 
A" computed for drip- line Cr isotopes with the {SLY4+REG- 
M} functional and plotted as a function of (e" — A"). Con- 
ventions for individual states are given in Fig. 1171 



section as it prevents any un-physical feature from ap- 
pearing with extreme surface pairing functionals. The 
properties of the last bound Cr isotopes are evaluated 
for different pairing functionals: (i) a = 1 and 77 S [0, 1], 
and (ii) rj — I and a £ [1,0.1]. 

Overall, neutron canonical energies evolve very little 
with (rj, a). The evolution of canonical pairing gaps A" 
with (77, a) is presented in Fig. [M] For surface-enhanced 
pairing functionals, well-bound orbitals residing in the 
center of the nucleus become less paired. On the other 
hand, pairing gaps of states close to the Fermi level in- 
crease as the effective pairing strength becomes more im- 
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Figure 22: (Color Online) Nhaio parameter for Cr isotopes cal- 
culated using renormalized REN-X pairing functionals with 
different form factors (see text). Results for regularized REG- 
X functionals diverge for very strong surface pairing correla- 
tions and are represented in dashed lines. 
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Figure 23: (Color Online) Same as Fig. [22]for 5_Rhaio. 



portant at the nuclear surface^''-'. Considering the the- 
oretical error bars, the values of iVhaio and (5i?haio pre- 
sented in Figs. [22] and [23] can be considered to be al- 
most independent of the density-dependent form factor 
of the pairing functional, although the anti-halo effect 
becomes more effective as a decreases. iVhaio and (5i?haio 
are maximal for the standard surface pairing functional 
•q = a — 1, for which the occupation of the 3si/2 state 
and its localization due to an increased coupling to the 
pairing field presents the most favorable compromise for 
the halo to develop. This reflects on the composition of 
the neutron halo in ^"Cr displayed in Fig. [23 Finally, 



* Keeping the neutron spectral gap unchanged in ^^"Sn, and as 
some orbitals become less paired, others display larger individual 
gaps. 



results obtained with the regularization scheme are also 
plotted in Figs. [22] and [23] for comparison. As long as 
a = 1, A^haio and (5i?haio are almost identical for the reg- 
ularization and renormalization schemes, which proves 
that results presented in Ref. Jj using the regularized 
scheme are valid and that both REG-X or REN-X func- 
tionals can be used safely with "standard" density de- 
pendencies. 
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Figure 24: (Color Online) Neutron single-particle pairing gaps 
A^ for ^''Cr and different REN-X pairing functionals. Con- 
ventions for labeling individual states are given in Fig. 1171 

As a conclusion, the impact of the low density char- 
acteristics of the pairing functional on halo properties 
is found to be small, as long as the adequate renormal- 
ization scheme is used to prevent the formation of the 
un-physical gas of bound di- neutrons. Consequently, ex- 
perimental constraints on pairing localization and the ef- 
fective pairing strength based solely on halo properties 
are unlikely. 



V. LARGE-SCALE PREDICTIONS 

We now turn to predictions of halo nuclei over the nu- 
clear chart. In addition, the dependence of such predic- 
tions on some of the characteristics of the Skyrme func- 
tional used in the calculations are briefly investigated. 



A. Results for different particle-hole functionals 

To study the impact of the particle-hole functional 
characteristics on the formation of halos, A^haio and 
<^^haio have been computed for Cr isotopes using the set 
of functionals presented in Sec. IIIIB II combined with a 
mixed-type REG-M pairing functional. The results ob- 
tained with SLy4 and presented in Ref. [ll] act as a ref- 
erence point. 

The halo parameters iVhaio and (5i?haio as well as neu- 
tron canonical energies are displayed in Fig. [2n]for ^"Cr. 
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Figure 25: (Color Online) Individual neutron contributions 
■/Vhaio.i for Cr and different renormalized REN-X pairing 
functionals. The contribution from the 2^5/2 orbital corre- 
sponds to vertically-dashed areas, and the one from the 3si/2 
state to filled ones. 



been used [39j . This particular series of parameteri- 
zations differ through their like-particle tensor cou- 
pling and an increasing spin-orbit coupling as one 
goes from T21 to T26. The like-particle coupling 
constant is (i) negative for T21, which corresponds 
to a repulsion between particles of identical isospin, 
(ii) zero for the functional T22, which makes it sim- 
ilar to SLy4, and (iii) increasingly positive from 
T23 to T26. All TXX functionals display other- 
wise the same infinite matter properties, as seen 
in Tab. IfVI the variation of single-particle prop- 
erties are solely due to the tensor and spin-orbit 
interactions. The results are shown in Fig. [26l As 
the like-particle tensor interaction becomes attrac- 
tive, both 3si/2 - 2(i5/2 and 2^5/2 - I59/2 energy 
gaps involving neutron states of the same parity 
decrease. The Fermi level is pushed up at the same 
time. Even though the core excitation energy scale 
E' is slightly decreased, the overall effect is domi- 
nated by the decrease of the scales E and AE which 
favors the formation of the halo [l|] . 



We observe that 

• Skyrme functionals with an isoscalar effective mass 
equal to one (m*l and T6) predict denser single- 
particle spectra around the Fermi level. As a con- 
sequence, the 2d5/2 — lff9/2 energy gap correspond- 
ing to the core excitation energy scale E' is reduced. 
Additionally, the Fermi level is shifted down, which 
enhances the separation energy E. Both effects 
contribute to hindering the halo formation. 

• The modification of K^o affects collective proper- 
ties, such as the breathing mode energy [9l|. ft also 
impacts the canonical spectrum. As a result, ^'^Cr 
is predicted to be unbound with the {SIII+REG- 
M} parametrization which has a particularly large 
incompressibility Kr^i — 355 MeV. 

• As the nuclear matter saturation density increases 
(from pI^^ to Psat)] fl^s nuclear interior becomes 
denser, as shown in Fig. [27l Through self- 
consistency, a denser nuclear interior generates 
a sharper surface, that ultimately makes weakly 
bound orbitals to be less coupled to the nuclear po- 
tential and thus even less bound. As a result, the 
density extends further out asymptotically. While 
the 2^5/2 — l<?9/2 energy gap remains the same, the 
tail excitation energy E decreases with psat and the 
halo factors increase. 

• Tensor couplings may also impact the formation 
of halos; as it was found for hght systems [Qgj . 
For medium-mass nuclei, a series of recent stud- 
ies have assessed the impact of tensor couplings on 
the evolution of sp herical single-particle shells [39, 
[Tool. [Toil Uol EqI. To study the impact on halos, 
newly available Skyrme T21-T26 functionals have 



By no means SLy4 or any other existing Skyrme 
parametrization of the nuclear EDF is to be seen as uni- 
versal. As exemplified above, it is found that the choice of 
the particle-hole Skyrme parametrization can affect sig- 
nificantly the neutron canonical spectrum close to the 
Fermi energy. As a result, the drip-line position and 
the halo energy scales are modified, hence the values of 
iVhaio and SRhaio- Scanning the set of parameterizations 
used in the present study, A'haio and (5i?haio can change 
by as much as 100%. Thus, solid experimental data on 
medium-mass halo nuclei will be useful in the far future 
to constrain some characteristics of the particle-hole part 
of the EDF. Note that this is at variance with the con- 
clusions drawn in Sec. IIV C] regarding the pairing part of 
the energy functional. 

Eventually, (reasonable) variations in the characteris- 
tics of currently used EDFs translate into large uncertain- 
ties in the prediction of halo properties. For instance, no 
halo nucleus would be predicted at the drip-line of Cr 
isotopes using the Skyrme SkM* parametrization, as the 
semi-magic nucleus ^''Cr is predicted to be the last bound 
Cr isotope in this case. The same care is to be consid- 
ered regarding the large-scale prediction of spherical halo 
candidates over the nuclear chart presented in the next 
section. 



B. Spherical medium and heavy mass nuclei 

Systematic predictions of halo properties from the 
{SLy4-|-REG-V} EDF are now presented. We restrict 
ourselves to even-even spherical nuclei, as predicted by 
the Gogny DIS interaction [104]. Among all even-even 
nuclei, we define the sub-set of "spherical" nuclei as those 
fulfilling the condition 
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Figure 26: (Color Online) Single-particle spectrum and halo properties for ^''Cr, computed with different particle-hole func- 
tionals (see text) and mixed-type pairing. Conventions for labeling individual states are given in Fig. 1171 



l/32N\ T7Z 



5 47r|Q2| 



167r iK^A 



<0.1 



(20) 



where Q2 is the axial mass quadrupole moment. The 
rather large interval allowed on (32 is arbitrary and does 
not distinguish between soft and rigid ground states. 
Such a condition provides a list of about 500 nuclei, in 
agreement with similar predictions made using Skyrme 
functionals [l05| , although the drip-line positions slightly 
differ in the two models. 

Results for A^haio are shown in Fig. [^51 We observe 
that (i) several isotopic chains are predicted to display 
neutrons halos, (ii) halos are predicted to exist only at the 
very limit of neutron stability, (iii) the maximum value 
of iVhaio is about ~ 0.7, (iv) very few heavy elements 
in the (Pt, Hg, Tl...) region are found to have a non- 
zero halo parameter iVhaio, (v) on the large scale, the 
halo phenomenon is very rare and almost accidental, (vi) 
looking at the best cases between Z = 2Q to Z = 100, 
the absolute and relative values of iVhaio decrease with 
nuclear mass. 

Results for (5i?haio are presented in Fig. [29] and confirm 
the above analysis on iVhaio- In particular, it is seen that 
the fraction of decorrelated nucleon has almost no influ- 
ence on the nuclear extension of massive nuclei. Only 



two very localized regions where the predicted halo sig- 
nificantly affects the neutron r.m.s. radius are found, e.g. 
for (i) Cr, Fe and Ni nuclei, and (ii) Pd and Ru isotopes. 
The drip-line isotopes of these elements are predicted as 
the best halo candidates for the {SLy4+REG-V} EDF. 

An analysis of single-particle properties of the best halo 
candidates is found in Fig. 1301 These nuclei have in com- 
mon the presence of very weakly bound s or p states. Al- 
though states with larger angular momentum contribute 
to the nuclear halo in some cases [l|, the presence of 
weakly bound ^ = 0, 1 orbitals remains mandatory for a 
significant halo to develop. That being said, no pure s- 
wave or p-wave halo has been found, which demonstrates 
the collectivity of the phenomenon in medium-mass sys- 
tems. 

The complementarity between the two criteria A^haio 
and (5i?haio appears more clearly through the large scale 
analysis. The plot presented in Fig.l^shows that the two 
observables are correlated within a given isotopic chains 
and the information carried by both quantities is some- 
what redundant in this case. On a larger scale however, 
the correlation pattern changes as the proton number in- 
creases (between Cr and Pd isotopes for instance), i.e. 
5R\iaio increases much less with A'haio with mass. 

Finally, one can turn to the particular case of Zr and 
Ca isotopes which have been predicted to be "giant halo" 
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Figure 27: (Color Online) Neutron densities for Cr 
(top: linear scale, bottom: logarithmic scale) with the 
{SLy4+REG-M} functional (reference), and refits of the 
Skyrme parameters corresponding to (i) a different isoscalar 
effective mass m* — 0.70 — » f.OO, and (ii) different nuclear 

0.145/0.160/0.175 fm"^ 



saturation densities: p^^ 



nuclei [H, \lM, [1Q3, EM \lM, EH- In the present 
study, the values of iVhaio and i5i?haio do not lead to 
such a conclusion. The first reason resides in the dif- 
ferent single-particle spectra predicted by the functional 
used in the pr esent study. Cons idering tha t relativistic 
models (TMl [uij, NLSH [HI, NL3 [lO^...) or non- 
relativistic one s (Sk I 4 Ill3ll...l predict weakly bound p 
states [H, Hoi, [lOl Il09l. Ill0l| at the neutron drip-line 
of those elements, it is likely that the application of our 
criteria on those results would lead to predicting the ex- 
istence of halos. However, we put into question the very 
notion of " giant halo" that comes from summing up the 
occupations of weakly bound orbitals. Such a counting 
procedure is qualitatively incorrect as nucleons occupying 
weakly bound orbitals are mostly located within the nu- 
clear volume. Such an unjustified counting is a reminis- 
cence of the denomination of "l(2)-neutron" halo used for 
light systems. However, the latter relates to the fact that 



these light nuclei are well described by a "core-|-l(2) neu- 
trons" cluster model whereas practitioners are well aware 
that only a fraction of nucleon resides in average beyond 
the classical turning point. We prove in the present work 
that the fraction of neutron that participate in the halo 
does not scale with A and remains in medium- mass nuclei 
of the same order as in light nuclei. The notion of spatial 
decorrelation is key to the meaningful definition of the 
halo region and the halo parameters iVhaio and (Si?haio- 



VI. CONCLUSIONS 

In Ref. [l|, a model- independent analysis method was 
developed to characterize halos in finite many-fermion 
systems. Eventually, the method leads to quantifying 
potential halos in terms of (i) the average number of 
fermions participating to it, (ii) its impact on the sys- 
tem extension. 

In the present paper, the versatility of the method is 
illustrated by applying it to several quantum mesoscopic 
systems, that is (i) few-body nuclei calculated through 
coupled-channels methods, (ii) atom-positron and ion- 
positronium complexes evaluated through the fixed-core 
stochastic variational method and (iii) medium-mass nu- 
clei calculated using the energy density functional (EDF) 
method. The ability of the analysis method to quantita- 
tively characterize halos of very different scales underlines 
the universal nature of the halo phenomenon, namely the 
quantum tunneling into the classically forbidden region. 

In the second part of the paper, the need for a coherent 
treatment of the ultra-violet divergence and of the den- 
sity dependence of the local pairing functional employed 
in EDF calculations is highlighted. Using a naive regu- 
larization scheme in combination with a strongly surface- 
enhanced pairing functional leads ultimately to the for- 
mation of a uniform gas of di-neutrons dripping out of 
nuclei. When the pairing functional is such that this spu- 
rious gas does not completely develop (most probably for 
numerical reasons) , it can be mistakenly interpreted as a 
halo. We point out that using a renormalization scheme 
or constraining the pairing functional through ab-initio 
calculations of pairing gaps in infinite nuclear matter pre- 
vents such problems from happening. 

In the third part of the present paper, the impact of 
pairing correlations on halo systems is further studied. 
First, it is shown that the pairing anti-halo effect might 
be counter-balanced by the pair scattering to less bound 
orbitals, possibly with small orbital angular momenta. 
The net effect of pairing on the formation of halos greatly 
depends on the detail of the single-particle spectrum of 
the nucleus under study. Second, low-i! orbitals are found 
to be less paired than high-^ ones, but not enough to 
really decorrelate from the pairing fieldj_as it is the case 
for (non-realistic) extreme conditions [42, 13 • Third, 
the low-density-dependence of the pairing functional is 
shown to have almost no effect on the formation and 
characteristics of halos, as long as it is combined with 
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Figure 28: (Color Online) A^haio parameter predicted by the {SLy4+REG-V} EDF for about 500 spherical nuclei. 
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Figure 29: (Color Online) 5_Rhaio parameter predicted by the {SLy4+REG-V} EDF for about 500 spherical nuclei. 



21 




^0.4 

^ 0.2 


0.0 


0.18 

o 

«0.12 

go.06 


0.00 








)K 



)K 



)K 



¥ ^ * 



)K 



^ 1 1 

A 


-i- 


~ 1 
1 1 


■ ! t ! 
1 1 


1 
1 
1 
1 
1 
T 


_ _ _, 


- : ^ ! 

A ® 1 


1 
1 
1 
1 


i ^ k 

* 5^ ® 



84 



Fe 



80, 



Cr 



88 



Ni 



8«Ni i^»Pd 



136 



Ru ''^Cr 
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{SLy4+REG-V} functional. The evolution of (5-Rhaio as a 
function of A^haio for Fe and Pd isotopes are highlighted. 



the renormalization scheme or constrained from ab-initio 
calculations, as mentioned above. 

On the contrary, halo properties significantly depend 
on the characteristics of the particle-hole part of the nu- 
clear EDF. Thus, medium-mass halo systems might be 
more useful to constrain the particle-hole part of the 
functional than its particle-particle counterpart. Unfor- 
tunately, such experimental data are not likely to become 



available any time soon, even with the next generation of 
radioactive beam facilities. Indeed, although the neutron 
density bec omes acce s sible thro i igh severa l experimental 
techniques [ui, [ul [ui, [u^, KM, [HI, extracting it 
in medium-mass drip-line nuclei is more than a challenge 
for the decade(s) to come. In addition, the neutron r.m.s. 
radius is not sufhcient to study halos quantitatively, and 
other probes have to be envisioned [6]. Of course, the 
precise determination of the neutron r.m.s. radius and 
associated neutron skin in non-halo systems is already 
crucial as i t provide s constraints on the physics of neu- 
tron stars 1 120 . Il2l| and on the nuclear symmetry en- 
ergy jl22l Il23j |. for instance. As a result, one should 
focus at first on the study of neutron skins in non-halo 
systems to constrain the isovector nature of the nuclear 
EDF. The fine tuning provided by extreme exotic sys- 
tems such as medium-mass halo nuclei will only come as 
a second step. 

Still, it is of theoretical interest to understand the 
structure of halo nuclei and to assess their occurrence 
over the nuclear chart. With that in mind, we per- 
formed large-scale calculations over all (predicted) spher- 
ical even-even nuclei. It was concluded that (i) several 
isotopic chains may display neutron halos, (ii) halos can 
only exist at the very limit of neutron stability, (iii) very 
few heavy elements, in the (Pt, Hg, Tl...) region, are 
found as possible halo candidates, (iv) on the large scale, 
the halo phenomenon is very rare and requires the pres- 
ence of a low-lying state with an orbital angular momen- 
tum £ < 1 at the drip line, (v) medium-mass halos are 
more collective than in light nuclei, with non-negligible 
contributions from several low-lying states, including or- 
bitals with i > 2, (vi) the impact of the halo on the 
nuclear extension decreases with the increasing mass. 

Let us now turn to potential works to come. First, 
several extensions of the newly proposed method can be 
envisioned. For instance, the extraction of halo proper- 
ties in deformed systems, which represent the majority of 
known or predicted nuclei, requires additional formal de- 
velopments, starting from a decomposition of the nuclear 
density in multipolar moments 



Pir] 



E 



pT^ir)Yrif). 



(21) 



The analysis regarding the relative asymptotic posi- 
tioning of spectroscopic amplitudes and the occurrence 
of crossing patterns should thus be adapted to multipo- 
lar moments /9™*(r). Effects such as the increased con- 
tamination of weakly-bound deformed orbitals by £ = 
components at the limit of stability, fav oring the forma- 
tion of halos in deformed systems [4J, 1 1241 . Il25l Il26| |. 
could be investigated. In addition to deformed systems, 
the method should be extended to odd-even and odd-odd 
nuclei. This requires to formulate the method for non- 
zero spin states. Also, the effects of an explicit treat- 
ment of long-range correlations, e.g. symmetry restora- 
tions and large-amplitude collective motion, on medium- 
mass halo nuclei should be studied in connection with 
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the analysis method proposed in the present work. The 
exphcit inclusion of such correlations would allow a more 
reliable extraction of A'^haio and SRhaio in hght systems 
through EDF calculations. For instance, the study of 
^^C, which is pre dicted to be a halo system by few-body 
calculations [l27| . could provide a bridge between few- 
and many-body techniques. 

Finally, it would be of interest to conduct the same 
study as the one presented in the present work using the 
Gogny effective interaction to further probe the depen- 
dence of the results on the characteristics of the energy 
density functional; e.g. finite-range effects. Such a study 
is not feasible using traditional codes making use of a 
harmonic oscillator basis because the asymptotic of the 
one-body density is not described correctly enough to 
extract A^haio and i5-Rhaio reliably. However, a newly de- 
veloped code expanding the quasiparticle states on the 
eigenstates of a Woods-Saxon potential should allow one 
to do so [12^. 
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